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Abstract 



. The factorization method is a convenient operator language formahsm for consider- 



ation of certain spectral problems. In the simplest differential operators realization it 
corresponds to the Darboux transformations technique for linear ODE of the second or- 
I der. In this particular case the method was developed by Schrodinger and became well 

. known to physicists due to the connections with quantum mechanics and supersym- 
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metry. In the theory of orthogonal polynomials its origins go back to the Christoffel's 
theory of kernel polynomials, etc. Special functions are defined in this formalism as 
the functions associated with similarity reductions of the factorization chains. 



^ • We consider in this lecture in detail the Schrodinger equation case and review some 

r~| , recent developments in this field. In particular, a class of self-similar potentials is de- 



scribed whose discrete spectrum consists of a finite number of geometric progressions. 
Such spectra are generated by particular polynomial quantum algebras which include 
g-analogues of the harmonic oscillator and su{l, 1) algebras. Coherent states of these 
potentials are described by differential-delay equations of the pantograph type. Appli- 
cations to infinite soliton systems, Ising chains, random matrices, and lattice Coulomb 
gases are briefly outlined. 
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1. Introduction 

From the very beginning quantum mechanics served as a rich source of good mathematical 
problems. It played a major role in the development of the theory of generalized functions, 
functional analysis, path integrals, to name a few. In the last two decades the "quantum" 
disease became so widespread in mathematics that it is difficult to guarantee that a randomly 
chosen mathematical term will not get its cousin with such an adjective in the foreseeable 
future (if it does not have already) . The influence of quantum mechanics upon the theory of 
special functions is also indispensable. A bright example is given by the angular momentum 
theory which has lead to the Racah polynomials. An aim of this lecture is to outline another 
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fruitful interplay between these two scientific fields inspired by the factorization method. 

This method was suggested by Schrodinger as a convenient operator language tool for 
working with quantum mechanical spectral problems [42] . It was reformulated as a problem 
of searching of the factorization chain solutions by Infeld [25]. The review [26] became a 
basic reference in this field. A fairly recent revival of the interest to this method occurred 
due to the discovery of its relation to the notion of supersymmetry (see, e.g. the review [22] 
and references therein). The author himself turned to this subject from the quantum chro- 
modynamics due to an idea of a generalization of the supersymmetric quantum mechanics 
[40] . Although it was clear that special functions play an important role in this formalism, it 
took some time to recognize that, heuristically, special functions are defined in this approach 
as the functions appearing from similarity reductions of the factorization chain (for a more 
precise formulation, see the Appendix). 

Prom the mathematical side, a related techniques for solving linear ordinary differential 
equations of the second order was proposed by Darboux long ago [12]. Its various gen- 
eralizations are referred to in the modern theory of completely integrable systems as the 
Laplace, Darboux, Backlund, dressing, etc transformations [1]; in the analysis of isomon- 
odromic deformations one deals with the Schlesinger transformations [29]. In the theory of 
special functions such transformations appear as contiguity relations. An important step 
in the development of the subject was performed in a series of papers by Burchnall and 
Chaundy [9], which contain even some parts of the operator formulation of the approach 
("the transference"). For a completeness, let us mention also the terms "shift operator" and 
"transmutation" , which are used in some other variations of the formalism. For a rigorous 
mathematical treatment of some aspects of the Darboux transformations technique or the 
factorization method, see [11, 14, 15, 30, 41]. 

As far as the discrete recurrence relations are concerned, actually, it is the Christoffel's 
theory of kernel polynomials that provides a first constructive approach of such kind to 
spectral problems. This theory is based upon the simplest discrete analogue of the Darboux 
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transformations. A complementary part to this Christoffel's transformation was found by 
Geronimus in [23]. First applications of the Schrodinger's factorization method to finite- 
difference equations is given in [37]. Recent developments in this direction are reviewed in 
[49]. Let us remark that the same techniques was rediscovered in the works on numerical 
calculations of matrix eigenvalues. More precisely, the well-known numerical LR, QR, g- 
algorithms, etc provide particular instances of the chains of discrete Darboux transforma- 
tions. The literature on the taken subject is enormous, it is not possible to describe all 
its branches. The list of references given at the end of this manuscript is not complete, it 
contains mainly the papers encountered by the author during his own work (additional lists 
can be found in [47, 49]). 

Despite of a vast variety of existing constructions, part of which was just mentioned, 
we limit ourselves in this lecture to the simplest possible case based upon the stationary 
one-dimensional Schrodinger equation 

Lip{x) = -ipxx{x) + u{x)ip{x) = Xip{x), (1) 

describing the motion of a non-relativistic particle on the line a; G R in the potential field 
u{x), which is assumed to be bounded from below. The operator L is called the Hamiltonian 
or the Schrodinger operator. For convenience, the particle's mass variable and the Planck's 
constant h are removed by rescalings of the coordinate x and the energy A. The equation (1) 
describes an eigenvalue problem for L and the eigenvalues (or the permitted bound states 
energies of the quantum particle) are determined from the condition that the modulus of the 
wave function ip{x) is square integrable, ip{x) e L^(R). 

Depending on the physical situation, the real line R may be replaced by an interval with 
an appropriate boundary conditions upon ip{x). However, if one is interested not in the 
spectra themselves, but in the differential operators L with some formal properties, then it 
is convenient to take x, A e C. 
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2. The factorization method 



Let us factorize the second order differential operator L — —(P/dx^ + u{x) as a product of 
two first order ones up to some real constant Aq: 

L = A+A- + Ao, ^^d/dx + f{x). (2) 

The function f{x) is a solution of the Riccati equation P{x)— fx{x)+\Q = u{x). Substitution 
of the ansatz f{x) — —(I)o,x{x)/(Pq{x) shows that —4>q,xx{x) + u{x)(f)o{x) — Xo(po{x), i.e. 4>o{x) 
is a solution of the original Schrodinger equation (1) for A = Aq. 

If f{x) is a smooth function, then is a hermitean conjugate of A' in -L^(R) and L 
is a self-adjoint operator. Under these circumstances Aq cannot be bigger than the smallest 
eigenvalue of L. Suppose that \q is the smallest eigenvalue of L, and let ipo{x) be the 
corresponding eigenf unction. It is well known that ipoi^) is nodeless and may be normalized 
to have the unit norm, ||'0o||^ — \ipo{x)\'^dx — L Then, 

/X fj^y 

where a, b are arbitrary constants, is the general solution of the equation L0o = Aq^q. For 
6 7^ 0, the resulting f{x) is singular at some point and the operators A"^ in (2) are not well 
defined. Let us exclude this situation, i.e. set a = 1, 6 = 0. 

Using the lowest eigenvalue eigenfunction of L, one can always factorize L as described in 
(2) with the well-defined operators A^. Vice versa, if one manages to find the factorization 
(2) such that the zero mode ■0o of the operator A", A't/jQ — 0, belongs to -^^(R) and f{x) 

is not singular, then Aq is the lowest eigenvalue and 

g- Jo f{y)dy 

= , p ^ 1/2 (3) 



is the corresponding normalized eigenfunction. 

Let us define now a new Schrodinger operator L by the permutation of the operator 
factors in (2) 

L = A-A+ + Ao, (4) 



whose potential has the form u{x) = P{x) + fx{x) + Aq. Evidently, one has the intertwining 
relations 

LA+ = A+L, A'L = LA-, (5) 

playing the key role in the formalism. Prom (5) one deduces that if iIj{x) satisfies (1), then 
the functions ijj = A~ip provide formal eigenfunctions of L. Indeed, 

Li)^L{A-%lj) = A-{L%Ij) = \{A-%1)). (6) 

Actually, this gives general solutions of the differential equation Li/i = Ai/i for all A, except 
of the point A = Aq, where the zero mode of A" is located. This problem is curable since 
one can find the general solution of the equation L0o = Aq^o separately: 

^°(^) = ^ + ^f 
where e are arbitrary constants and -00 is the lowest eigenvalue eigenfunction of L (note 
that A+i^Q^ = 0). 

Denote as A„, A„ and ipn, i^n — ■^~'4^n discrete eigenvalues and corresponding eigenfunc- 
tions of L and L respectively. From (6) it follows that the spectra A„ and A„ almost coincide 
An = A„, n = 1, 2, . . .. The only difference that may occur in these spectra concerns the zero 
mode of the operator A~ . In fact, the point A = Aq does not belong to the spectrum of L. 
Indeed, since ■00 £ -^^(R-), it follows that iJjq'^ ^ -^^(R), and, as a result, for any g,e the 
function (7) cannot be normalizable. 

Thus, Ai is the lowest eigenvalue of L and the point Aq was "deleted" from the spectrum 
of L. Repeating the same procedure once more, i.e. taking L = A'^A~ + Ai and permuting 
the operator factors, one can delete the point A = Ai, etc. This procedure allows one to 
remove an arbitrary number of smallest eigenvalues of L. Often it is much easier to find 
the smallest eigenvalue of a given operator than the other ones. If the lowest eigenvalue of 
L is determined separately by some means, then, by construction, it will coincide with the 
second eigenvalue of L, etc. This observation is the central one in the factorization method 
[42] since it reduces the problem of finding complete discrete spectrum of a taken operator 
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to the problem of finding lowest eigenvalues of a sequence of operators built from L by the 
"factorize and permute" algorithm. 

One can invert the procedure of a deletion of the smallest eigenvalue. Namely, a given L 
with known lowest eigenvalue Ai may be factorized as L = A~A'^ + Aq, with Aq < Ai. If the 
zero mode of A~ is normalizable, then the operator L — A'^A" + Aq has the same spectrum 
as L with an additional inserted eigenvalue at an arbitrary point A = Aq. If one factorizes 
L (or L) in such a way that f{x) is a non-singular function, but none of the zero modes of 
the operators A"^ are normalizable, then the discrete spectra of L and L coincide completely 
(an isospectral situation). 

There are more complicated possibilities for changing spectral data of a given Schrodinger 
operator L. E.g., if the first factorization is "bad" , in the sense that L has a singular potential 
and considerations given above are not valid, then one may demand that after a number of 
additional refactorizations one gets a well defined self-adjoint Schrodinger operator. In this 
way one can delete not only the lowest eigenvalues of L or to insert the new ones, but delete 
or insert a bunch of spectral points above the smallest one [30]. In particular, two step 
refactorization procedure allows one to imbed an eigenvalue into the continuous spectrum, 
etc [14]. 

Let us give now a "discrete time" formulation of the construction. Denote L = Lj and 
L = Lj+i and take j G Z {j may be treated as a continuous parameter and we could denote 
L = Lj_i — all this is a matter of agreement). This gives an infinite sequence of Schrodinger 

operators Lj = —(P/dx^ + Uj{x) with formal factorizations 

Lj = A+AJ + \j , Af^^ d/dx + fj (x) . (8) 
Neighboring Lj are connected to each other via the abstract factorization chain 

Lj+, = A+^,AJ^, + A,+i = AjA+ + A,-. (9) 
Intertwining relations take the form 

Lj — Lj^iAj , LjA^ — Aj'Lj+i. 
7 



Substituting explicit forms of Af into (9) one gets a differential-difference equation upon 

ifjix) + + ff{x) - ff^,{x) = i,j ^ A,.+i - A,.. (10) 

This chain was derived in [25] and the problem of searching "exactly solvable" spectral 
problems was formulated as a problem of the search of solutions of (10) such that the points 
Xj define the discrete spectrum of an operator L, say, L = Lq. E.g., one may try to find 
solutions of the equation (10) fj{x) in the form of power series in j {Xj are considered as 
unknown functions of j). As shown in [25, 26] the finite term expansion occurs iff fj{x) — 
a{x)j + b{x) + c{x)/j, where a,b,c are some elementary functions of x. This leads to the 
2F1 hypergeometric function and well known "old" exactly solvable potentials of quantum 
mechanics. 

As evident from the construction, the constants Xj,j > 0, determine the smallest eigen- 
values of Lj under the condition that the zero modes of Aj are normalizable and fj{x) 
are not singular. In general, the j ^ j + 1 transitions may describe all three possibilities 
— removal or insertion of an eigenvalue and isospectral transformations (sometimes it is 
convenient to parameterize inserted eigenvalues as Aq > Ai > . . . > A„). 

For a positive integer n, let us introduce the operators 

Mj = Aj+n-l • • • ^j+l^j ' = ^t^t+l ' ' ' 

The intertwining relations 

guarantee that for almost all A solutions of the equations Lj-^^^ — Xt/j^^^ and = 
Xi/jU+n) g^j^g related to each other as V^^^"^ <^ Mfi/j^^ and i/j^^ oc M/V^^'^"^- As a result, 
the product M~ should commute with Lj and Mj M'j' should commute with Lj^n- A 
simple computation leads to the equalities 

n— 1 n—1 

Ml Mr =\[{L,- A,+,), Mr Ml = [] ( - A,+.). (11) 

k=0 k=0 



Let ^^^\x) e L2(R) be an eigenfunction of Lj with the eigenvalue A and some finite norm 
If we set = M^ip^^\ then 



||V^(^-+")|p = (A-A,)---(A-A,+._OI|V^(^-) 



Zeroes on the r.h.s. for A = A^ for some k indicate that the corresponding eigenvalues were 
deleted from the spectrum. Under the condition that zero modes of Aj,j = 0, . . . ,n — 1, 
are normalizable and nodeless, AJi/jq^ — 0, ||'0o''^|| — 1? oiie finds that the functions 



define the unit norm eigenfunctions of Lq with the eigenvalues A„: Lot/j^^ — Xnfpn^ and 



The main advantage of the factorization method consists in its pure operator language 
formulation. One may replace the Schrodinger operator by any other (higher order differen- 
tial, finite-difference, integral, etc) operator L admitting factorizations into the well defined 
operator factors of a simpler nature A"^ (for some applications they are not necessarily her- 
mitean conjugates of each other). In all cases one deals with the abstract operator factoriza- 
tion chain (9) with an appropriate interpretation of the constants Xj. It should be noted that 
this method does not have straightforward generalizations to the multidimensional spectral 
problems, only some of its features are preservable, see e.g. [3] . 

3. Supersymmetry 

Supersymmetry is a symmetry between bosonic and fermionic degrees of freedom of par- 
ticular physical systems. The corresponding symmetry algebras are distinguished from the 
standard Lie algebras by the presence of anticommutator relations. The simplest superalge- 
bra is realized upon two neighboring operators in the factorization chain. Define the matrix 




l^f 11 = 1- 
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operators (supercharges) 



A+ 




Q- 




where are the factorization operators defined in (2). They form the following algebra 
(see, e.g. [22]) 

{Q+, Q-} = H, [H, g±] = (g±)2 = 0, (12) 



where the Hamiltonian 

/ 



H 



L-Xo d 



2 



y L-Xo 



describes a particle with spin 1/2 (a fermionic variable) on the line in an external magnetic 
field, (T3 is the Pauli matrix. In terms of the hermitean supercharges Qi = + , Qi = 
{Q~^ — Q~)/i, the algebra takes the form {Qi, Qj} = 2HSij, [H, Qj] = 0. 

As a formal consequence of the algebra (12) the discrete spectrum of H is doubly degen- 
erate (such degeneracies are characteristic to the supersymmetry) with possible exception of 
the lowest eigenvalue which cannot be negative (this is another constraint upon the super- 
symmetric systems). 

This simple construction was generalized in [40] to a symmetry between particles with 
parastatistics. The corresponding symmetry algebras are polynomial in the generators. E.g., 
in the simplest case involving a parafermion of the second order (or the spin 1 particle) H 
is given by a 3 x 3 diagonal matrix with the entries Lq, Li, L2, the corresponding symmetry 
generators being {Q^)ij = Af6ij-i and its conjugate. For the hermitean charges (5i,2 defined 
as above, one has now cubic relations 

QiiiQj, Qk} - '^H6jk) + cyclic perm, of k = 0, [H, Qi] = 0. (13) 

The discrete spectrum of H is now triply degenerate with possible exception of two smallest 
eigenvalues. 
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One can go further and propose other modifications of the algebra (12). For instance, 
one can g-deform it [46]: 

g+g- + q^Q-Q+ = H, (g±)2 = o, hq^ = q^^q^h. 

This hfts the degeneracy of spectra, but creates a nontrivial scahng relation between spectral 
points. As a result, it opens a way for building various g-harmonic oscillator models. Another 
construction proposed in [4] refers to a polynomials generalization of (12). It appears after 
a reduction of the algebra (13) and similar higher order polynomial relations to a two- 
dimensional subspace of eigenf unctions: 



{g+, g-} = p^{H), [H, g±] = {Q^f = o. 



where Pn{H) = UkZliH - Xk) and 
/ 

g+ = 






/ 




\ 


H = 






J 


\ 





Ln J 



g-- 



This is a pure supersymmetry again in the sense of a symmetry between bosons and fermions. 
However, the consequences are quite different from the standard case (12). In particular, 
there are no such severe constraints upon the smallest eigenvalue of H (the vacuum energy) . 

All these general constructions are useful because they contain new concepts providing 
nonstandard viewpoints upon physical systems and new mathematical tools for their explo- 
ration. Special functions emerge when one starts to work with particular "exactly solvable" 
models with such symmetries. 



4. Darboux transformations 

Let us describe the Darboux transformations technique in its modern appearance. Consider 
compatibility conditions of the following three linear equations 

Lji)^^\x, t) = \ij^^\x, t), Lj = -dl + Uj{x, t), (14) 
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= Ajilj^^\x,t), AJ = + fj{x,t), (15) 

(x, t) = Bji/j^^^ {x, t) , Bj = -Adl + 6uj{x, t)d^ + ?>Uj^^{x, t) , (16) 

where t is some additional continuous parameter ("evolution time") and Uj{x, t), fj{x, t) are 
some free functions. The compatibihty condition of (14) and (15) generates the intertwining 
relation AjLj — Lj+iAj , the resolution of which yields the constraints Uj — /J — fj^x + ^j, 
Uj+i — Uj + 2fjx, where is an integration constant. Thus one arrives again to the infinite 
chain of nonlinear differential-difference equations (10) and the formal factorizations (8). 

The compatibihty condition of the equations (14) and (16) yields the celebrated Korteweg- 
de Vries (KdV) equation: 

Uj,t{x, t) - 6uj{x, t)uj^x{x, t) + Uj^xxxix, t) = 0, (17) 

a completely integrable Hamiltonian system with an infinite number of degrees of freedom 
[1]. If we express Uj{x,t) through fj{x,t) in the equation (17), then it takes the form 

{-dx + 2fj) Vj{x, t) = 0, Vj{x, t) = fj^t - Q{f- + Xj)fj,x + hxxx- 

It appears that the compatibility conditions of this equation with the Darboux transforma- 
tion (15) imposes an additional non-trivial constraint Vj = 0, which is called the modified 
KdV equation. 

A substitution of fj — —4>ol/(po^ into the relation uj — /| — fj^x + yields — ^o^lx + 
Uj(f)Q^ — Xj(f)Q^ and the Darboux transformation (15) takes the form: 

v^(^-+^) = (rfVi^-)-0gv'W)/#, 

where (pQ^ is a particular solution of the j'-th equation in the sequence (14) for X = Xj. 

Let Lj(j)ll^ = A,-+fc4^'\ i.e. let 0?(a;,t) are formal eigenfunctions of Lj with the eigenvalues 
Aj+fc equal to constants of integration mentioned above for all k. Denote as W{(f)i, . . . , 4>n) — 
det(9^"~^0fc) the Wronskian of a set of functions 0^. Using the identities VF(^(x)0i, ^(x)02) = 
eix)W{(f)i,(f>2) and 
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one can show that [11] 

uj+n{x,t) = u,{x,t) - 2a2logW^(#, . . . ,rf2i), (18) 

wUi'\...,(PiiA 

fM^-d.lo,^^^^^^, (19) 

i.^^-^--^'\x, t) = ) ° ;: • • ' " ; ^ = {a^ + /,+„) ...{d^+ f,)^p^^\x, t). (20) 

These formulae give an explicit representation of the transformed potentials Ujj^n in terms 
of the initial Schrodinger equation solutions 0^^^. 

Since the parameter i is a dummy variable during the Darboux transformations, in fact 
one deals simultaneously with different solutions of the KdV equation: if Uj{x,t) satisfies 
(17), the same is true for Uj+n{x,t). This gives a way for building new complicated explicit 
solutions of the KdV equation starting from the simple ones. For instance, one may take 
Uo{x, t) — which gives (l)^\x, t) — a^e'^'^^-^'^^^^ + fc^e"'*'^^'*"^'^"^*, and for some special choice 
of the signs of hm/chm the potential Un{x,t) becomes a nonsingular refiectionless potential 
with n discrete spectrum points = —t^ln < 0- It defines the famous n-soliton solution of 
the KdV equation. 

In the formulae given above one may permute any pair of the (f)^^ functions and the final 
result is not changed, i.e. such a permutation of the Darboux transformations is a symmetry 
of the factorization chain (10). In the simplest case, it has the following explicit form [2]: 

Jk — Jk 7 r~F' •'^+1 ~ ■'^+1 ~'~ ~1 ^ ~ '^k+l-i '^k+l — -^k, 

Jk+l + Jk Jk+1 + Jk 

all other fj{x), Xj staying intact for j ^ k, k+1. This discrete symmetry may be "discretized" 
further by passing to the discrete Schrodinger equation (or the three term recurrence relation 

for orthogonal polynomials) and the corresponding discrete Darboux transformations [48]. 

The well known KdV tau-function is introduced as Uj{x,t) = — 20'^ log Tj{x,t). In its 
terms relations (18), (19) are rewritten as Tj+„ = W{(f)Q \ . . . , 4>n-i)'Tj, fj — —dxlogrj+i/rj. 
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Introducing the variables pj = —dx logrj, one can write Uj = 2pj^x and fj = pj+i — pj. As a 
result, the relation between uj and fj yields the equation 

(Pj+i + Pj)x - {Pj+1 - Pjf = Aj-, 

which starts to play the role of the factorization chain. The function tj is a very convenient 
object since its zeros in x correspond to poles of the potential. 

5. Operator self-similarity and quantum algebras 

Let us turn now to the problem of searching particular solutions of the operator factorization 
chain (9). At first glance it is not clear how to proceed, but the harmonic oscillator problem 
— a base model for the whole quantum mechanics and quantum field theory — provides a 
guiding idea. One has to try to form from the factorization operators Af some nontrivial 
symmetry algebras. In fact, the relations (11) look already as defining relations of some 
algebra, but they are not closed — the relations between the eigenvalue problems for Lj 
and Lj+n are too weak (they are valid for any starting potential Uj{x)). In order to close 
the system, one has to assume that there is an additional relation between Lj and Lj+n, 
say Lj+n — g{Lj), which would force the operators Mf" to map eigenfunctions of a taken 
operator Lj onto themselves. 

In the simplest case one demands that the sequence of Hamiltonians Lj is periodic: 
Lj_|_7v = Lj for some period > 0. As a result, commute with Lj, [M^,Lj] = 0. 
This is a remarkable fact since the existence of additional conserved quantities may simplify 
solution of the eigenvalue problem for Lj. In the differential operator realization of Af 
this leads to the commuting differential operators [9] . There is a generalization of this pure 
periodicity condition to the periodicity up to a twist condition Lj+jv = ULjU~^, where U 
is some invertible operator. This leads again to commuting operators [Bf,Lj] = 0, where 
Bf = Mj~U, Bj = U~^M^ , but now there is an essential additional freedom in the choice 
of U. 
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Another possible "closure" or a reduction of the sequence of operators Lj consists in the 
requirement of their periodicity up to a constant shift and a twist, Lj^]\f = ULjU~^ + /i, 
where is a constant. This results in the ladder relation [Lj, B^] — i/i-Bj*^ and the operator 
identities 

Af-l N-l 

Bj'Bj = Y[ {Lj - \j+k), BJB^ = Jl (Lj + /X - Xj+k), 

k=0 k=0 

where B^ operators are defined as in the previous case. Denoting B"^ = Bq,L = Lq, one 
can form a polynomial algebra 

[L, B^] = ±fiB^, [B+, B-] = Pn-i{L), (21) 

where Pn-i{x) is a polynomial of the degree N — 1 in x. For a representation theory of such 
algebras, see e.g. [45]. Note that for = 1 this is the Heisenberg-Weyl or the harmonic 
oscillator algebra, and for iV = 2 it coincides with the su{l, 1) algebra. 

Quantum algebras, or gr-analogues of the algebras (21) appear from the following operator 
self-similarity constraint imposed upon the chain (9): 

Lj+N = q^ULjU-^ + /X. (22) 

When the numerical factor q'^ ^ 1, one can remove /i by the uniform shift Lj Lj+ji/ (l—q^). 
Therefore we assume below that = 0. Substitute (22) with fj, = into (11). Then the 

operators L = Lq, B^ = MqU, B~ = U^^Mq satisfy the following identities 

N-l N-l 

LB^^q^'^B^L, B+B- = ll{L-Xj), B' B+ = H {q"" L - Xj) . (23) 

j=0 j=0 

For N — 1 these relations provide a reahzation of the g-harmonic oscillator algebra 

B-B+ - q^B+B- = p, [B^, p] = 0, (24) 

with p — Ao(?^ — 1) . One can set p = 1 by taking the normalization condition Aq = l/{q^ — l). 
Such a g-analogue of the Heisenberg-Weyl algebra was encountered in physics long ago 

[10, 18]. In the modern times it became quite popular due to the inspirations coming from 
the quantum groups, see e.g. [35]. For N = 2 relations (23) determine a particular g-analogue 
of the conformal algebra su{l, 1) admitting the Hopf algebra structure, etc [46, 47]. 
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6. Self-similar potentials 



Consider the Schrodinger equation realizations of the algebraic relations described in the 
previous section. Let us start from the closure L = ULU~^ + fj,, where L = —d^ + u{x) 
is the Schrodinger operator in the notations of Sect. 2 and t/ is a translation operator 
Uf{x) — f{x + a). As a result, the operators take the form S+ = {—dx + f{x))U, B~ — 
U~^{dx + f{x)) and one gets the Heisenberg-Weyl algebra [L, B^] — -^jiB^, [B~ , B'^\ — ji. 
The potential entering L is defined as u{x) = p{x) — fx{x) + Aq, where f{x) satisfies the 
following nonlinear differential-delay equation 



For a — one gets from this equation the standard harmonic oscillator model f{x) oc 
x,u{x) oc x'^, which is related to the Hermite polynomials. Although for a ^ and /i = 
the author has found a meromorphic solution of (25) in terms of the Weierstrass P-function 



in general it is quite difficult to build its solutions analytic in some region, especially for 
11 ^ 0. Let us remark, that the equation (25) may be derived also after imposing the 
constraint fj+i{x) — fj{x + a) upon the chain (10). 

Turn now directly to the relations (23). In this case one can take U as the dilation (or g- 
difference) operator: Uf{x) — \q\^^'^f{qx). For real q ^0 this is a unitary operator W — U"^. 



Taking L = Lq, = U-\dx + /jv-i) ■ ■ ■ (5. + /o), B+ = (-9, + /o) ■ ■ ■ {-dx + /^-i)?7 



one realizes the identities (23), provided fj{x) satisfy the following system of nonlinear 
differential-delay equations: 



{fix) + f{x + a)), + f{x) - f{x + a)=ii. 



(25) 



1 V'{x - xq) - V'{a) 

2 V{x - xo) - r{a) ' 



{h{x) + h{x)). + /o (^) - = /^o, 



{fN-i{x) + qfo{qx))x + fN-i{x) - q^ fl{qx) = Hn-i- 



(26) 
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Equivalently, these equations appear from the chain (10) after imposing the following con- 
straints 



having a simple group-theoretical interpretation. 

One easily arrives at (27) using the Lie's symmetry reduction technique (for a differen- 
tial context of this theory, see e.g. [38], and for an extension to difference equations, see 
e.g. [31, 36]). First, one notices that if fj{x), jij are solutions of the chain (10), then their 
discrete scaling transformation qfj{qx), q^/ij gives solutions of (10) as well. Analogously, 
a shift in the numeration fj{x) fj^N{x), fJ^j /J'j+N maps solutions of (10) into so- 
lutions. Let us consider a set of self-similar solutions of (10) which is invariant under a 
combination of these two symmetries. E.g., demanding that these two transformations are 
equivalent to one another, we arrive to (27). For A^" = 1, this reduction may be rewritten as 
fj{x) — fo{q-'x), \j = q^^Xo in which form it was first met in [43]. The general g'-periodic 
reduction (27) was found by the author [46] in an attempt to g-deform the parastatistical 
supersymmetry algebras (13). 

Solutions of the equations (26) have a quite complicated structure and general methods 
of solving such differential-delay equations give relatively weak results. E.g., for A" = 1 in 
[44] the existence and uniqueness of solutions analytical near the x — point was proved 
under some constraints, and in [32] existence of the nonsingular for ,x G R solutions was 
demonstrated (the ^ oo asymptotics of such solutions is not determined yet completely). 
Some understanding of the complexity of general solutions is obtained from considerations 
of special values of the parameters when fj{x) can be expressed through some known 
functions. 

Let q be arbitrary and ff{x) = Y^mZl^m - J2m=o l^m- This yields Lq = -d'^/dx^ or 
the free quantum mechanical particle, which acquires in this way a g-algebraic interpretation 



Suppose that fo{x) is not singular at x — 0, then in the crystal base limit — > the 





(27) 



[47]. 
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potential uq{x) boils down to the general KdV A^-soliton potential. 

Substitute now q = I into (26) and assume that E^=o/^m 7^ 0- Then for N = 1,2 one 
easily gets the potentials Uo{x) oc x^,ax^ ^-h/x^. For N — 3 the corresponding system of 
equations for fj{x) provides a "cyclic" representation of the Painleve-IV equation [8, 50]. 
The N — A case is related to the Painleve- V function [2] , etc. This gives also a commutator 
representation of the Painleve equations which was noticed first in [16]. For q = —1 one 
gets similar situations with the functions fj{x) obeying certain parity symmetry [47]. If 
Z]m=oA*m — ^^^^ ^ some cases of even N the potential Uo{x) is expressed 

through the hyperelliptic functions [50, 51]. 

The case when g is a primitive root of unity, q'^ — l,q ^ ±1, is quite interesting [44]. 
For any odd n and, under certain conditions, for even n potentials are expressed through 
hyperelliptic functions with additional crystallographic or quasi-crystallographic symmetries. 
For example, for N — 1 and q^ — 1 or q"^ — 1 one recovers the Lame equation with the equi- 
anharmonic or lemniscatic Weierstrass P-f unctions. Thus, this classical differential equation 
of the second order appears to be related to the representations of the g-harmonic oscillator 
algebra for q a root of unity. 

The general family of self-similar potentials unifies some of the Painleve functions with 
hyperelliptic functions. Due to the connections with quantum algebras, it may be taken as a 
new class of nonlinear g'-special functions (viz. continuous gr- Painleve transcendents) defined 
upon the differential equations. The standard basic hyper geometric functions [21] appear in 
this formalism through the coherent states. 

As far as the discrete spectrum of self-similar potentials is concerned, from the factoriza- 
tion method and (27) it follows that in the simplest case it is composed from N independent 
geometric progressions: XpN+k — \q^^, k — 0,...,N — 1, peN. The same follows from 
the theory of unitary representations of the algebra (23). The only condition for the validity 
of this formal conclusion is that the functions fj{x) are non-singular for a; G R and positive 
for X — > -|-oo and negative for x — > — oo. 
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7. Coherent states 



There are many definitions of quantum mechanical coherent states [39]. In the context of 
spectrum generating algebras, such as (21) and (23), they are defined as eigenfunctions of the 
corresponding lowering operators. Such eigenfunctions play the role of generating functions 
for the space elements of irreducible representations of a taken algebra. Let us describe 
briefly coherent states for the quantum algebra (23) (for more details, see [47]). 

Let us denote as |A) eigenstates of an abstract operator L entering (23), L|A) = A|A). 
Suppose that the operators are conjugated to each other with respect to some inner 
product (c|A). We assume that the spectrum of L is not degenerate. Then the action of 
operators upon |A) has the form: 

JV-l N-l 

i^-|A) = n iX-Xk\Xq-'), B^\X) = n VV-A,|Ag2). 

k=0 k=0 

Let N be odd, < < 1 and Aq < ... < Aat-i < g^Ao < 0. Then for A < the 
operator L may have a discrete spectrum formed only from up to N geometric progressions 
corresponding to the lowest weight unitary irreducible representations of the algebra (23). 
This follows from the fact that B~ is the lowering operator for the A < states and Y\k=d{'^~ 
Afc) becomes negative for A < Aq. Since B~\Xk) — 0, this problem is avoided for special 
values of A, namely, for A = Xp^^k = X^q"^^, p G N. For normalizable \Xk), one gets series of 
normalizable eigenstates of the form |Apjv+fe) oc {B^Y \ Xk). 

Coherent states of the first type are defined as eigenstates of 

B-\a)^!:^ = a\a)^!l\ k = 0,...,N-l, (28) 

where o; is a complex variable. Representing |q;)?^ as a sum of the states |Apjv+fe) with some 
coefficients, one finds 

oo / D D \ 

p=0 V Ao/ Afc, Aat-i/ Afc / 

where j^ipj^_i is a standard basic hypergeometric series with the operator argument z = 
(— 1)-^q;S+/Ao • • • Ajv-i and C^^\a) is a normalization constant (in the set of parameters 
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of the N^N-i function the term \k/\k = 1 is assumed to be absent). The superscript k 
simply counts the number of lowest weight irreducible representations of the algebra (23) 
each of which has its own coherent state (or a generating function). The states \(y)^- are 
normalizable for jap < |Ao • • • Aiv-i|- Coherent states of this type were first constructed for 
the gr-harmonic oscillator algebra (24) in [5], for a particular explicit model of them, see e.g. 
[6]. 

It is not difficult to see, that zero modes of the operator L define a special degenerate 
representation of the algebra (23). Since B"^ operators commute with L in the subspace of 
these zero modes, one can consider them as coherent states as well. 

Unusual coherent states appear from the non-highest weight representations of (23) cor- 
responding to the A > eigenstates of L. Let A+ > be a discrete spectrum point of L. 
Then the operators generate from |A_|_) a part of the discrete spectrum of L in the form 
of a bilateral geometric progression A+g'^",n e Z. Since nfe^o^(''^+ ~ "^k) > for arbitrary 
N and A+ > 0, the number of such irreducible representations in the spectrum of L is not 
limited. Moreover, a continuous direct sum of such representations may form a continuous 
spectrum of L. 

An important fact is that for A > the role of lowering operator is played by B~^. 
Therefore one can define coherent states as eigenstates of this operator as well [47] : 

B+\a)+^ a\a)+. (29) 

Suppose that the bilateral sequence A+g^" > belongs to the discrete spectrum of L. Then 
the states \a)+ are expanded in the series of eigenstates |A+g^"'): 

n=-oo \Ao/A+, . . . , AiV-l/A+ / 

where Qip^ is a bilateral g-hypergeometric series with the operator argument z — aB~ / (— A+)^ 
and C{a) is a normalization constant. In this case the states \a)+ are normalizable if 

IcKp > |Ao ■ ■ ■ A;v— 1|- 
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Suppose that the A > region is occupied by the continuous spectrum. The correspond- 
ing states I A) may be normaUzed as (cr|A) = XS{X — a). Let A'" = 1 and the states |A) are 
not degenerate. Taking p = 1 in (24), we get the expansion 

ds)_r^^^^ r A'^^|A)dA 



\a) 



C{a) f 
Jo 



where 

_ 27ris - \R{aq'^y/l - q^) 

Is — '■, n ) S & L. 

m q-^ 

There is a countable family of such coherent states which have a unit norm for |q;|^ > 
1/(1 — q^) under the following choice of the normalization constant C{q): 

' ^ ^' io (-Ag2(l-g2).^2)^ 3in^^ (^2.^2)^ 

where v — \n. \aq\/l — q^\/ Ing. The integral under consideration is a particular subcase of a 
Ramanujan gr-beta integral [21]. 

We were considering until now the abstract coherent states. In the case of Schrodinger 
equations with the self-similar potentials one has the doubly degenerate continuous spectrum 
for A > 0. The structure of \<y)± states in this case is very complicated. Even for the zero 
potential case (the free particle) we get a highly nontrivial situation. Let L — —d'^/dx^ and 



B- = U-\d/dx+ 1/^1 -q^), B+ = {-d/dx+ 1/^1 -q^)U, 

where U is the dilation operator, Uf{x) = q^^'^f{qx), < g < 1. It can be checked directly 
that B-B+ - q^B+B- = 1 and L = B+B' - 1/(1 - q^). The equation B-'il^-{x) = a'il)-{x) 
coincides now with the retarded pantograph equation, which was investigated in detail in 
[28] . As follows from the results of this paper, in the retarded case the pantograph equation 
does not admit solutions belonging to L^(R). However, the operator S+ has infinitely 
many normalizable eigenstates. The equation B^iIj^^^x) = aip'^{x) coincides now with the 
advanced pantograph equation 

^^^{x) = -aq-'/'^^{q-'x) + -^i=^^+(a;), (30) 
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which has infinitely many solutions belonging to iv^(R) [28]. For the free particle realization 
of the N > 1 symmetry algebras, one gets the generalized pantograph equations considered 
in [27]. As one of the open problems in this field, let us mention a need to characterize a 
minimal set of solutions of the equation (30) providing a complete basis of the Hilbert space 



8. Solitons, Ising chains, and random matrices 

Consider the KdV A^-soliton solution UN{x,t) — — 29^ logT]v(a;, i), where tn is a Wronskian 
of N different solutions of the free Schrodinger equation (18). Actually, there are many 
determinantal representation for this tau-function, e.g. [1] 



The variables ki are known to describe amplitudes of sohtons (they are related to the eigen- 
values of Ln as Aj = — /c|/4), Of^ /ki are the t — Q phases of solitons and kf are their 
velocities. 

A subclass of the self-similar potentials appearing from the g-periodic closure (27) can 
be viewed as a particular infinite sohton potential. Indeed, consider the KdV pM-soliton 
solution with the parameters kj subject to the constraint /cj+M = ?%,0 < q < 1, and take 
the p ^ oo limit. The limiting potential has the self-similar spectrum kpM+m = Q^km,fn = 

0. . . . , M — 1, and an infinite number of free parameters 6f\ 

The scaled potential u{x, t) — q'^u{qx, qH) has the same solitonic interpretation with 
the phases 9j(x,t) — 6j(qx,qH) — kj+MX — k^jj^j^t + Of^- Let us demand that ^j+M — 

1. e. impose the additional constraint 9j{qx,qH) = 6j+M{x,t), which is seen to match with 
the condition (27). In this picture the discrete dilation x qx,t q^t just deletes M 
solitons corresponding to the smallest eigenvalues of L — this feature is the simplest physical 
characterization of the self-similar potentials. 





(31) 



i,j = l,2,...,Ar. 



22 



Let us turn now to the statistical mechanics appUcations. A^-sohton tau-function of the 
KdV equation (31) can be represented in the following Hirota form [1, 24]: 

/ N-l \ 

tn^ ^ exp\ AjCTiCTj + ■ (32) 

<Ti=0,l \0<i<j<iV-l i=0 J 

The coefficients Aij, describing phase shifts of the i-th and j-th solitons after their collision, 
have the form 

(/Cj + fcj)^ 

As noticed in [33] , for 9^ — this tau-function coincides with the grand partition function 
of a lattice gas model on a line (for a two-dimensional Coulomb gas picture, see [34] and the 
next section). In this interpretation the discrete variables ai describe filling factors of the 
lattice sites by molecules, 9'^^^ is a chemical potential, and are the interaction energies of 
molecules [7]. 

It is well known that a simple change of variables in the lattice gas partition function 
leads to the partition function of some Ising chains: 

N-l 

Zn^ Y e"'''', E= Y J^JS^SJ - E HiSi, (34) 

Si=±l 0<i<j<N-l i=0 

where N is the number of spins Sj = ±1, J^j are the exchange constants, is an external 
magnetic field, P — 1/kT is the inverse temperature. Indeed, a replacement of filling factors 
in (32) by the spin variables Ui — (sj -|- l)/2 yields 

1 1 

^ i<j ^ 3=0 



where 



1 . 1 1 



/5J,, = --yl,„ f3H, = -9, + - Y Ar (36) 

Similar relations with Ising chains are valid for the whole infinite Kadomtsev-Petviashvili 
(KP) hierarchy of equations and some other differential and finite difference nonlinear in- 
tegrable evolution equations. Corresponding tau-functions have the same form (32) with 
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a more complicated structure of the phase shifts Aij and phases 6i (a partial list of such 
equations can be found in [1, 24]). 

The general relation between A'"-soliton solutions of integrable hierarchies and partition 
functions of the lattice statistical mechanics models established in [33] brings a new point of 
view upon the self-similar potentials as well. Namely, as shown in [33] self-similar spectra 
can be deduced from the condition of translational invariance of the infinite spin chains. 
E.g., let us demand that the spin chain is invariant with respect to the shift of the lattice 
by one site j ^ j + i which assumes that Jj+i^+i = . As a result, the exchange constants 
Jij (or the soliton phase shifts A^j) depend only on the distance between the sites \i — j\- 
This very simple and quite natural physical requirement forces ki to form one geometric 
progression (N.B.: uniquely) 

ki = koq\ q = e~^°, Aij = 2 In | tanh Q;(i - j) I , (37) 

where ko and < g < 1 are some free parameters. In a more complicated case one demands 
the translational invariance with respect to the shifts by M lattice sites, i.e. Ji+M,j+M — Jij, 
and this results in the general self-similar spectra kj+M = Qkj. The main drawback of these 
relations between the solitons and Ising chains is that for fixed q the temperature T (or (3) 
is fixed too, which is clearly seen from the comparison of (36) and (37) (in the Coulomb gas 
interpretation one actually finds that (3 — 2). 

For finite chains < j < — 1, the translational invariance is not exact. The infinite 
soliton limit N ^ oo corresponds to the thermodynamical limit in statistical mechanics. In 
this picture, all coordinates of integrable hierarchies {x,t and higher hierarchy "times") are 
interpreted as particular parameters of the external magnetic field Hi. Since < g < 1, 
the X and f-dependent part of Hi decays exponentially fast for i — > oo. As a result, in the 
limit A^ — > oo only the constants Of*^ are relevant for the partition function (speaking more 
precisely, they determine the leading asymptotics term of the partition function at A^ — oo). 
Although in the thermodynamic limit the x, t-dependence is washed out and the Ising chain 
is periodic, presence of an infinite number of free parameters Of*^ does not allow one to find 
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a closed form expression for the leading term of Z^. However, for the self-similar infinite 
soliton systems, characterized by the periodicity of Of'^ or, equivalently, by the periodicity 
of the external magnetic field H^^m — H^, such expressions do exist. 

Take M — 1 and let the magnetic field be homogeneous Hi — H. Then in the KdV 
equation case one gets an antiferromagnetic Ising chain: < |tanha;(i — j)| < 1 and 
Jij = —Aij/4f3 > (a similar picture is valid for the general M). The exchange has a 
long distance character but its intensity falls off exponentially fast and, as a result, phase 
transitions are absent for nonzero temperatures. However, for a — > or — > 1 the exchange 
constants are diverging. Under appropriate renormalizations one can build a model with 
a very small non-local exchange such that a phase transition takes place even for non-zero 
T. 

In order to compute the partition function one can use determinantal representations for 
the tau- function, e.g. the Wronskian formula (18). As shown in [33], for the translationally 
invariant Ising chains in a homogeneous magnetic field tau-functions become determinants 
of some Toephtz matrices and their natural M-periodic generalizations. E.g., for M = 1 the 
following result has been derived: Zj^ — >■ exp(— for N — >■ oo, where the free energy 
per site // has the form 

-/3fj{q,H) = \n^^-Lh_^+ d,,\n{\p{,.)\'-qtanh'(3H), (38) 

(g2; q^)oL 47r Jo 



(g4e-;g4)2^(g4e-...^4)2^ 4sin2(z//2) ^el{u/2,q^) 
(^1,4 are the standard Jacobi theta-functions). In the derivation of this formula the Ramanu- 
jan iipi sum was used at the boundary of convergence of the corresponding infinite bilateral 
basic hypergeometric series. A standard physicists' trick was used for dealing with that, 
namely, a small auxiliary parameter e was introduced into the original expression for the 
density function p(i/) 

. . ^ {q';q')l e^-'-' ^ (gV^-; g^)oo(g^e-+-; g^).o 
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which guarantees the absolute convergence of the series for e > 0. Ahhough the hmit e — > 
leads to some singularity in v on the right hand side of (39), it does not lead to divergences 
after taking the integral in (38) (a more precise mathematical justification of this physical 
"harmlessness" is desirable). The total magnetization of the lattice m(if) = —dnfi = lim 

N^oo 

J2^jQ^{si) has the following appeahng form 

In order to change the fixed value of the temperature, one may try to replace (37) by 
Aij = 2n In I {ki — kj) / {ki + kj) | , where n is some sequence of integers, and look for integrable 
equations admitting A'"-soliton solutions with such phase shifts. This is a highly nontrivial 
task, and in [33] only one more permitted value of the temperature was found. It corresponds 
to n = 2 and appears from special reduction of the A'"-soliton solution of the KP-equation 
of B-type (BKP) [13]. The general BKP r-function generates much more complicated Ising 
chains than in the KdV case. The corresponding exchange constants have the form 

^ ^ (a, + a,)(6, + 6,)(a, + 6,)(6, + a,)' ^^'^ 

where a^, hi are some free parameters. Translational invariance of the spin lojtticG, J^ij — 
J{i — j), results in the following spectral self-similarity 

ai = q\ k = bq\ g = e'^", (42) 

where we set ao = 1 and assume that a > 0. This gives the exchange 

tanh^a(^-j)- (6 -1)7(6 + 1)^ 
^ 4 coth^ a{i - j) -{b- 1)7(6 + 1)2 ■ 

The parameter b is restricted to the unit disk |6| < 1 due to the inversion symmetry b — > b"^. 
For —l<b<—q one gets now the ferromagnetic Ising chain, i.e. Jjj < 0; for g < 6 < 1 
and b = e"^ ^ —1 the chain is antiferromagnetic, i.e. Jij > 0. In the thermodynamic 
limit N ^ oo the partition function Zjv is given again by the determinant of a Toeplitz 
matrix, which is diagonalized by the discrete Fourier transformation. Similar picture holds 
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for arbitrary M-periodic chains, see [33] for explicit expressions for the magnetization and 
technical details of computations. 

It is well known that the n-particle lattice gas partition functions define probability 
distribution functions in the theory of random matrices. In a natural way, partition functions 
of the Ising models considered in the previous section may be related to the probability 
distribution functions of some random matrix models with a discretized set of eigenvalues. 
Let us describe briefly this correspondence comparing the KP A'"-soliton solution and the 
Dyson's circular ensemble. 

Dyson has introduced an ensemble of unitary random nxn matrices with the eigenvalues 
€j = e^'^\ j — 1, . . . ,n, such that after the integration over the auxihary "angular variables" 
the probability distribution of the phases < 0j < 27r takes the form P(i0i . . . d(f)n oc 
Y{i<j |ej — ej\'^d(j)i . . . d(j)n. One may relax some of the conditions used by Dyson and work 
with a more general set of distribution functions. In [19] Gaudin has proposed a circular 
ensemble with the probability distribution law 

2 

#1 . . . #„, (43) 



Pd0i . . . d0„ oc J\ 

i<j 



^3 



containing a free continuous parameter cu. It interpolates between the Dyson's case {cu — 0) 
and the uniform distribution (a; = 1). Let us note in passing, that the paper [19] contains 
in it implicitly a specific model of the g-harmonic oscillator. The model (43) admits also an 
interpretation as a lattice gas on the circle with the partition function 

/•27r r2TT f \ 

Z^GC ... o?0i . . . d0„ exp V V{(f)i - (pj) , 
Jo Jo \ I 

where (3—1/ kT = 2 is fixed and the potential energy has the form 

/3m-*.) = ln(l + -.^|^;^), ^ = e-^. (44) 

A surprising fact is that the grand partition function of this model may be obtained in 
a special infinite soliton limit of the A'^-soliton tau-function of the KP hierarchy [33] . This 
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means that the finite KP sohton solutions provide a discretization of the model, namely, 
they define a lattice gas on the circle. This leads to random matrices with a discrete set of 
eigenvalues. E.g., one may take unitary n xn matrices with eigenvalues equal to A'"-th roots 
of unity, i.e. ej — ex^{2mmj/N), rrij — 0, . . . , N — 1. The probability measure is taken to 
be continuous in the auxiliary "angular" variables of the unitary matrices and discrete in 
the eigenvalue phase variables (pj. More precisely, the integrals over (pj are replaced by finite 
sums over rrij and the continuous model is recovered for rrij, N ^ oo with finite rrij/N: 

n N-1 iv-i „27r 



^ ' mi=0 m„=0 



The n-particle partition function becomes 

n N-1 



r2'K /'27r 
/ d(pi... d(j)n 

Jo Jo 



(45) 



27rY 

Iv) 



N-l 

EE n 

m„=0 l<i<j<n 



mi=0 



€i 6 



while the grand canonical ensemble partition function takes the form 



Z{U!, I 



N-l 



n=0 



ni 



E E ^mkCTmCrk + {0 + V) J2 



(46) 



=0,1 



\0<m<A;<Ar-l 



m=0 



where rj — ln(27r/A^) is an excessive chemical potential, and 

sin^(7r(m-A:)/A^) ^ (a^ - ak){bm - h) 
sin^(7r(m - k)/N) + sinh^ 7 (a^ + hk){hm + Ofc) ' 

are the KP sohtons phase shifts for a special choice of the parameters 



In 



—ujam-i m = 0, 1, . . . , A?" — 1. 



(47) 



To conclude, the grand partition function of the Gaudin's circular ensemble coincides 
with the particular infinite soliton KP tau-function at zero hierarchy "times". The root of 
unity discretization of circular ensembles has been considered by Gaudin himself [20], where 
a connection with Ising chains was noticed as well, but the relation with soliton solutions 
of integrable equations was not established. The BKP hierarchy of equations suggests a 
generalization of the distribution law (43) to 



Pd0i . . . d0„ oc Jl 



60 ~t~ €j 



Si + uej 



d(j)i... dcpn, 



(48) 
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and its discrete analogue, which were not investigated yet appropriately. 

It is well known that the n-tuple Selberg integral provides an explicit evaluation of the 
n particle Coulomb gas partition functions for arbitrary values of the inverse temperature (3 
(i.e. not just for /3 = 1, 2, 4, corresponding to orthogonal, unitary and symplectic ensembles). 
It would be interesting to know whether a similar universal formula exists for the grand 
partition functions of lattice Coulomb gases. 

9. Lattice Coulomb gas on the plane 

Solution of the Poisson equation on the plane, AV{r,r') — — 27r5(r — r'), defines the elec- 
trostatic potential V{z, z') = — In |^ — ^'|, created by a charged particle placed at the point 
z' — x' + iy' . In the bounded domains with dielectric or conducting walls, the potential 
has a more complicated form since the normal component of the electric field E = — VV^ 
should vanish on the surface of a dielectric, = 0, while the tangent component is zero 
at the metallic boundary, Et — 0. For simple geometric configurations of boundaries an 
introduction of artificial images of charges may simplify solution of the Poisson equation. 

The energy of an electrostatic system ( "plasma" ) consisting of particles in a bounded 
domain of the plane is 

l<i<j<N l<i<N l<i<N 

where zj = Xj + iyj and qj are the particles' coordinates and charges. The first term is the 
standard Coulomb energy, the second one describes an interaction with the boundaries (or the 
charge- image interaction), and the last term corresponds to the external fields contribution. 

Suppose that our plasma is composed from the particles of equal charges, qj — +1 (in the 
two component case qj — ±1), upon a discrete lattice F on the plane. The grand partition 
function of such lattice Coulomb gas is 

n=o "'■ zier z„6r 
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where jj, is an effective chemical potential. In general case one can rewrite Gjv in the form 
[7] 

G^- ^ exp(i^W^(^,^')<^(^M^') + E^(^M^))' (50) 

W {z, z') = -Pq{z)q{z')V{z, z'), w{z) = f,{z) - (5 {q\z)v{z) + q{z)(t>{z)) , 

where a{z) = or 1 depending on whether the site with the coordinate z is empty or 
occupied by a particle, and the functions q{z) and iJ,{z) characterize distribution of particles of 
different types. E.g., in the two-component case, when q{z±) — ±1 charges occupy the {z±} 
sublattices, one has iJi{z±) — ii±. Note that these 2D lattice gases describe simultaneously 
some 2D Ising magnets. 

Let us take now the A^-soliton tau-function (32) and replace in it the soliton number j 
by a variable z taking N discrete values. As a result, it may be rewritten in the form 

riv= E exp(i^A,,,a(zM/) + ^^^(z)a(^)). (51) 

<t(z)=0,1 ^ z^z' {z} 

But this is precisely (50) — one just needs to identify the phase shifts A-^z' with the Coulomb 
interaction potential W(yZ, z'), and the phases 9{z) with the function w{z). As a result of this 
observation made in [34], one can construct a number of new solvable models of Coulomb 
gases in addition to already known ones (see, e.g. [17, 19, 20] and references therein). 
For the KP-hierarchy sohton solutions one has [13] 

where tp is the p-th KP hierarchy "time" and a^, are some arbitrary functions of z. For the 
BKP-hierarchy A,,> are given by (41) and e{z) = e^^\z) + E^i(a^^"^ + bf-^)t2p-i, where 
t2p-i are the BKP evolution "times". 

Let us take — z — x + iy, — —z* — —x + iy. Then in the KP case 

A^^,^W{z,z') = -2V{z,z') =21n|^-^'| -21n|2*-^'|, (53) 
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where V{z, z') is the potential created by a positive unit charge placed at z' over a conductor 
with its surface occupying the y = line. In this case V{z, z') solves the Poisson equation 
with the boundary condition St{y = 0) = 0. The same potential is created by a positive 
charge placed at the point z' and its image of opposite charge located at the point {z')*. 

Similar to the random matrix models situation, the correspondence between 2D lattice 
Coulomb gases with specific boundary conditions and solitons is valid only for fixed j3, which 
is found from the comparison of (53) with (50) to he (3 = 2. Since w{z) = 0{z), one finds an 
exphcit expression for the zero time sohton phases 9^^\z) in (52): 

e'^°\z) = fi- (3{ln\z* - z\ + (l){z)), /3 = 2, (54) 

where the middle term corresponds to the charge-image interaction, and (f){z) is the potential 
created by a neutralizing background of some density p(r), A0(r) = — 27rp(r), (p^iy = 0) = 0. 
The harmonic term Z^^ol-^^^l-^*)^)^? — —I3(t>ext{z) in (52) corresponds to an external electric 
field. One may conclude that the imaginary hierarchy times evolution of the KP soliton 
solutions describes the behavior of a 2D lattice Coulomb gas in a varying external electric 
field. 

Usefulness of the conformal transformations z f{z) is well known in the potential 
theory. With their help one can map plasma particles to various regions. So, the map 
z — > puts them inside a corner with the vr/n angle between the conducting walls. The 

exponential map z eJ^^I^ leads to a gas inside the strip ^ <'^z < L between two parallel 
conductors, etc. 

As to the BKP equation case, the choice — z,bz — z* leads to the Coulomb gas placed 
inside the upper right quarter of the plane with the dielectric and conductor walls along the x 
and y-axes respectively. The discrete temperature appears to be the same, /3 — 2. A curious 

set of solvable dipole gas models is generated by the self-similarity constraints imposed upon 
the spectral data of the corresponding multi-soliton systems. For further details, see [34]. 
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10. Appendix. A heuristic guide to special functions 
of one variable 



There are many handbooks and textbooks on special functions. However, none of them 
contains a formal list of properties which a function should have in order to be "special". 
One usually talks on functions of some particular type (hypergeometric, L- functions, etc). 
R. Askey has suggested one universal definition that if a function is so useful that it starts to 
bear some name, then it is "special" . Another essential, but not so universal, characteristic 
refers to the asymptotic behavior. Namely, for special functions one is expected to be able 
to deduce asymptotics at infinity from a known local behavior, i.e. the connection problems 
should be solvable — such an approach is advocated by people working on the Painleve type 
functions. 

To the author's taste these two definitions are relying upon the secondary features of 
special functions. One has to have already the functions in hands in order to start to 
investigate their properties. If one takes as a goal the search of new special functions, then 
it is necessary to find a definition containing a more extended set of technical tools for work. 
In this respect it should be stressed that even the term "classical special functions" appears 
to be not so stable. E.g., it is by now accepted that the family of classical orthogonal 
polynomials consists of not just the Jacobi polynomials and their descendants, but of the 
essentially more complicated Askey- Wilson polynomials invented just two decades ago. 

The group theory provides a number of tools for building new functions, but, unfor- 
tunately, connections with the representation theory often provide only interpretations for 
functions already defined by other means. Still, the symmetry approach seems to be the 
central one in the theory of special functions. In particular, the key "old" special func- 
tions appear from separations of variables in the very simple (and, so, useful and universal) 
equations [38]. For the last decade the author's research was tied to the following working 
definition: special functions are the functions associated with self-similar reductions of spec- 
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tral transformation chains for linear spectral problems. Speaking differently, special functions 
are connected with fixed points of various continuous and discrete symmetry transformations 
for a taken class of eigenvalue problem. This definition works well only for special functions 
of one independent variable (which may, however, depend on an infinite number of param- 
eters) and even for them it does not pretend to cover all possible special functions. On the 
one hand, this definition comes from the theory of completely integrable systems, where a 
search of self-similar solutions of nonlinear evolution equations is a standard problem [1]. 
On the other hand, particular examples found from this approach show that it has in its 
heart the contiguous relations — linear or nonlinear equations connecting special functions 
at different values of their parameters. 

Schematically, a heuristic algorithm of looking for these "spectral" special functions con- 
sists of the following steps: 

1. Take a linear eigenvalue problem determined by differential or difference equations. 

2. Consider another linear equation involving variables entering the first equation in a 
different way and having the same space of solutions. 

3. Resolve compatibility conditions of the taken linear problems and derive nonlinear 
equations for free functional coefficients entering these problems. When the second equa- 
tion is differential, one gets the continuous fiows associated with the KdV, KP, Toda, etc 
equations. When the second equation is discrete, one gets a sequence of Darboux transforma- 
tions performing some discrete changes in the spectral data and providing some discrete-time 
Toda, Volterra, etc chains. 

4. Analyze discrete and continuous symmetries of the latter equations in the Lie group- 
theoretical sense, i.e. look for nontrivial continuous and discrete transformations mapping 
solutions of these nonlinear equations into other solutions. 

5. Consider 5e//-simi/ar solutions of the derived nonlinear equations, which are invariant 
with respect to taken symmetry transformations. As a result, some finite sets of nonlinear 
differential, differential-difference, difference-difference, etc equations are emerging which de- 
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fine "nonlinear" special functions. Solutions of the original linear equations with coefficients 
determined by these self-similar functions define "linear" special functions. 

The last two steps are heuristic since, despite of a big progress in the general theory 
of symmetry reductions (see, e.g. [31, 36]), no completely regular method has been built 
yet. E.g., the reductions used in the derivation of the associated Askey- Wilson polynomials' 
recurrence coefficients [49] and in the discovery of new explicit systems of biorthogonal 
rational functions described in Zhedanov's lecture at this meeting did not find yet a purely 
group-theoretical setting. 

Another important ingredient of the theory of special functions, which was not listed 
in this scheme, is the theory of transcendency. Painleve functions are known to be tran- 
scendental over the differential fields built from a finite tower of Picard-Vessiot extensions 
of the field of rational functions. In solving differential (difference, or whatever) equations 
one has to determine eventually which differential field a taken solution belongs to (e.g., 
whether it belongs to the field of functions over which the differential equation is defined). 
As an open problem in this field, which was discussed partially in [47], we mention a need 
to find a differential (or difference) Galois theory interpretation of the self-similar reductions 
of factorization chains. 
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